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Polymer's network is treated as an anisotropic fractal with fractional dimensionality D = 1 + e 
close to one. Percolation model on such a fractal is studied. Using real space renormalization group 
approach of Migdal and Kadanoff we find threshold value and all the critical exponents in the perco- 
lation model to be strongly nonanalytic functions of e, e.g. the critical exponent of the conductivity 
was obtained to be e~'^exp(— 1 — 1/e). The main part of the finite size conductivities distribution 
function at the threshold was found to be universal if expressed in terms of the fluctuating variable 
which is proportional to a large power of the conductivity, but with e-dependent low-conductivity 
cut-off. Its reduced central momenta are of the order of e'^^"^ up to the very high order. 

Disordered structures with fractional dimensionality arise in many various physical applications B. In particular, 
a number of papers are devoted to the processes of diffusion on fractals (|-^ . In our resent paper Q a new type of 
fractals is introduced: nearly one- dimensional strongly anisotropic ones. There we were dealing with the problems of: 
i) percolation and ii) variable range hopping on them. One motivation to study these problems were the experimental 
data on the structure of some classes of oriented conducting polymers, and on their unusual conducting and dielectric 
properties |^-|^. The other is their conceptual significance, and, in particular, the possibility to deal directly with the 
distribution functions (DF) of strongly fluctuating random variables. We extend here an approach, developed in |^ 
in such a way, that we are able to obtain the DF of finite sample conductivity on the percolation threshold. 

In Q conducting polymers were modelled as fractal oriented networks. Namely, we assume all chains to be ori- 
ented in one spatial direction, and coupled transversely through various size perfectly conducting islands. Fractional 
dimensionality was defined as follows: In an L-size cube chains form a set bundles, connected within this cube. If 
the cross-section of maximal bundle scales as L*^ for large enough L, where < e < 2, then we have d* ~ 1 + e 
-dimensional network. Obviously e = for purely one-dimensional systems (sets of disconnected chains) . The main 
feature of the fractals, constructed from oriented Id chains is their self-similarity: the system at any scale looks like 
subdivided into bundles, which are in turn subdivided into smaller ones, etc. Our hypothesis here is that oriented 
polymers network structures are of this type (with d* = 1 + e close to 1, e ^ 1), at least in some wide enough length 
scales interval, e.g. from the scale of polymer's fibrils (hundreds of nm) down to molecular scales. Transmitting 
electronic micrographs (see e.g. Q]) seems to confirm this hypothesis. 

A regular example of such a fractal is the hierarchical lattice [^O[0 , constructed by the infinite repetition of two 
steps: i) connection of n bonds in sequence to form n-bond chain, and: ii) assembling of m n-chains in parallel to 
form m-bundle, which is treated as a new bond on the next stage, etc. The fractal dimensionality of such a system 
is: d* = him/ Inn -\- 1. An exact real space renormalization group may be written for hierarchical lattice, which 
becomes the renormalization group of Migdal and Kadanoff (RGMK) ||lO|-|l|l if to set n ^ 1, m ^ 1 with the value 
of d* fixed. Of course, the real polymers structures are not regular ones, and the requirement of self-similarity here 
is to be treated in statistical sense. Nevertheless, we shall use the RGMK scheme. The additional argument for 
using this approach in nearly Id case is that the RGMK is exact in one dimension, therefore one may hope to obtain 
meaningful results when the dimensionality is close to 1. This method was applied to the percolation conductivity 
problem two decades ago by Scott Kirkpatrick He had found values of critical exponents of correlation length 
and conductivity near threshold, using the RGMK equations for conductivity momenta near threshold, truncated at 
the first moment. Though he had not considered explicitly the case of dimensionality close to one, this method, if 
appropriately applied, gives the right dependence of conductivity exponent on e up to pre-exponential factor. 

The RGMK method may be formulated in a quite simple phenomenological fashion. Suppose we have some random 
d-dimcnsional medium with the local conductivity — fluctuating random variable. Let us consider the A-size cube 
within the medium. Its conductivity Sa, and resistivity R\ — l/S^, are fluctuating random variables with some DFs, 
defined (in the Laplace representation) as: 

F(a,A) = (exp(-aI](A))); Q(s, A) = (exp (-si? (A))) . (1) 

If we change the size of the cube, A — + A' = tiA, we arrive at some new random variables S(A'), R{X'), the DFs 
are also changed, of course. The cube's enhancement may be treated as a combination of n-times expansions in 
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one "longitudal" (arbitrarily chosen) spatial direction, and in cZ — 1 "transverse" ones. Further we shall transit to 
infinitesimal transformation, so the order of operations is not important. For finite length rescaling we assume that, 
changing the size n times in the longitudal direction, we arrive at n resistivities in sequence, with the resulting 
(specific) resistivity: 
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where Ri are assumed to be independent random variables. Therefore, we have: 

Q{s) = ^cxp (-si?)) = Q" (s/n) , 

and for infinitesimal transformation with n = 1 + 6X/X, the variation of the DF is: 
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where Si means the variation due to longitudal rescaling. Quite similarly, after the transverse rescaling, when the 
cross-section is enlarged m = n'^~^ times, and the resulting conductivity is assumed to be the arithmetic average of 
m statistically independent ones, we have the conductivities DF P{a) transforming in the same manner as in (H), but 
with n replaced by m. For infinitesimal transformation m = 1 + e6X/X, where e = d — 1, and we have for transverse 
rescaling: 



StP{<J,X) = 
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Using the integral identity: 
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where Ji is the Bessel's function, one can write the relation between conductivities and resistivities DFs in the form 
of Hankel's transformation: 



Q(s, A) = 1 - x/i / — Ji(2V^)F 
Jo 



(s,A) 
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The reverse transformation is quite the same. Now we are able to write down both transverse and longitudal variations 
in terms of either conductivities or resistivities DF. Adding both variations for e.g. conductivity DF we arrive at the 
following evolution equation upon size rescaling: 
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This equation should be completed with Eq. (y) to form a closed set. 

We may introduce the probabilities of A-cube to be disconnected (i.e. to have zero conductivity, or infinite resistivity) 
c(A) which, taking into account the definitions of DFs dll), may be written as: 



c(A) = P (+00, A) = 1 - g (+0, A) 



Assuming in Eq. (^ cr = +00, we have: 
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This equation has three fixed points: two stable ones, c — and c = 1, corresponding to connected and disconnected 
systems resp. in the thermodynamic limit, and the unstable fixed point, c — ct, < ct < I, 
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ect In ct = (1 - ct) In (1 - ct) , 
corresponding to the percolation threshold. The correlation length exponent is given by: 



dc 



e + 1 + ln Ct + ln(l - a) 



For nearly-ld systems, e ^ 1, we have: 

Cf = i/=l/e. 
It is possible to rewrite Eq.(^) using WKB-type approximation, assuming: 

P(a) =c+(l-c)exph(/.(a)] , 
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with 0(0) = 0, and (ct) ±oo as ct ^ ±oo rapidly enough (more rapidly then ±-y/|CT[, as we shall see later). The 
important point also is to assume analyticity of P (ct) and of Q{s) at least within some finite width stripe along the 
real axis. Using the relations: 
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where ' are the Hankel's function of first and second kind, resp., and directing the cut of the function 
h[^\2^/z) / ^ in z complex plane along positive real half-axis we may replace the integrals with Ji-function along 
the positive real half-axis in Eqs. (||) and (^) with the ones containing h[^\ along the following contour C: from 
-l-oo — iO to S — iO along the bottom shore of the cut, then from 6 ~ iO to 6 + iO along the almost closed anticlockwise 
(5-circle, and finally from 6 + iO to +oo + iO along the top shore of the cut. Thus we have: 
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Assuming |s| to be large enough, one may replace li\ in the latter integral by its asymptotic expression 
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and to treat this integral in the saddle point approximation. Afterwards, the same procedure may be performed with 
the integral in Eq.(|6|). As a result, we have the evolution equation in the saddle point or "WKB" approximation to 
be: 
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It seems that saddle point approximation is valid only at large enough values of ct. The other approximation for 
the evolution equation is possible if to set: P (ct) = c + {1 — c) e""' -I- V' (ct), and to linearize Eq.(^ with respect to ^. 
The remarkable fact is that after substitution of the above expression into Eq.(p^), we arrive at the same linearized 
equation for ij) (ct, A). Thus, we have some reason to look for an appropriate solution of Eq.(p^) in the whole complex 
plane ct. 

At the percolation threshold, c = ct, one may to look for the solution of the RG evolution equation in the form: 
P(ct, A) = P(ctA^°), where a — t/v, t and ly are critical exponents of the conductivity, S cx (q — c)*, and of correlation 
length, ^ cx |c — Ct| . Then the equation ( |l^ ) becomes an ordinary differential one of the second order. It appears 
to be more convenient to use the function (p{x) = — In [(P(a;) — q) /(I — Ct)] instead of P{x). Introducing (/)o = 
— In [cf /(I — Ct)], we have: 
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= [l + e - a) xcf)' - (f) + e[g {(j) - (po) - g (-00 ) 



(15) 
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where (/(c/i) = (e^ + l) In (l + e*) . An equation for 0o which follows from Eq. was used in the derivation of Eq. 
(|l5|). The latter may be easily solved after the substitution: 

z{<j)) = exp[-2{l + e-a)x(j)'] , (16) 

with (j) — new independent variable. Requiring z{(j)) as (/> faster then exp [— (1 + e — a) 0], we have: 

/•oo 

z(0) = (l + e-a)e-(i+^-")* / dy exp{- {1 + e - a) y + 2e[g {y - cl>o) ~ g ■ (17) 

Jtf, 

The normalization condition 0(0) = implies 2;(0) = 1, from which it follows that: 

/>oo 

(1 + e-a)/ dyexp{-{l + e-a)y + 2e[g{y-cf>o)-g{-cf>o)]}^l, (18) 



which is an equation for a. 

Comparing the values of a, obtained by the solution of Eq.(^8|), and by the numerical investigation of the evolution 
of the evolution equation (||) Q one can see, that both methods give the same results at any dimensionality. This, 
together with the considerations presented above, prompts us to consider the saddle point solution as an exact one. 
Of course, the RGMK method itself is an approximate one. In e.g. three dimension we have from Eq. |l^: a « 1.891. 
On the other hand, the best possible at present numerical results ||l^ give a — 2.25 ± 0.04. So, the RGMK method 
may be not very bad even for 3-d systems. 

In case e ^ 1 our previous result ||] : 

a= — exp ( , t = a/£, (19) 



may be easily reproduced from Eq. ( [l^ ) . In case of e ^ 1 one may obtain from Eq. ([T 



a = e--e"^ (20) 



Finally, the function (/)(a;) may be determined as a reverse of the equation 

'1 + e-a 1 



Cx = (pcxp 
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The arbitrary integration constant C corresponds to arbitrary choice of the unit of conductivity, or, alternatively, of 
the length scale at the threshold point. 

Thus the conductivities DF in the normal representation, n(S(I], A) = Ct + (1 — Ct)II [(A/Aq)" S], where the scaling 
hmction II(y) may be expressed as: 

n(y) = / exp [xy - <l)ix)] = - / exp [-</. + yx(</))] , (22) 

the last equality was obtained through the integration by parts. However, there is some trouble when evaluating 
integral in Eq. (p2[). Namely: the function x{(j)) is singular at = ^„ = (j)Q + i'K{2n+l), n- integer. The origin of these 
singularities is in the procedure of analytic continuation made when the RGMK approach was formulated. It can be 
illustrated as follows: Let us assume the initial distribution of conductivities to be: Poia) = c+{l—c)e~"' .After putting 
m identically distributed conductivities in parallel, the Laplace of the DF of their sum, Pi{a) = [c+ (1 — c)e~'^]"^ , 
has TO-th order zeroes at an = (Tq + i7r(2n + 1), n-integer, cto = ln[(l — c)/c], which turns into singularities after 
analytic continuation to nonintcger to. This points us, that it is the procedure of the transition from integer rescaling 
factor transformation (which is exact for a hierarchical structure) to the infinitesimal one (which no explicit structure 
corresponds to) that is the reason of these singularities. So, these singularities are to be treated as artiflicial ones, 
and should be avoided during the integration in Eq. (p^). In general, this restricts our knowledge of the DF 11 (y) 
with its low- and large-conductivities asymptotic behavior. 

At large scaled conductivities y, shifting integration contour in Eq. ( p^ ) to the region Sft^ > 00, one has the following 
asymptotic expression for the DF: 
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Shifting the integration contour in Eq.(|22|) to the region > (/ig, we arrive at the foUowing expression for the n(?;) 
DF at small y-region: 
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with: 
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More detailed results are available in the limiting case e <C 1. Here we have at 3fi0 < — l/fi; the following 
expression for x{(j)), truncated at the first order of Ct and of a: 



lnx(0) = ln0 + Cf 



°> - 1 - ^ 



Evaluating Taylor series of <f) {x) at a; = 0, we obtain central momenta of the conductivity to be of the order of a 



(27) 



a + Cf 
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On the other hand, using in Eq.(p2|) the asymptotics of x[(f)) at 5ft0 < 0, ^ 1, which may be justified at 
sufficiently large y, we have, after the proper change of the integration variable: 



(29) 



where 02 = ^/(l + ^Y' ~ c-t, and the new fluctuating variable was introduced: 

f7 = e«^^y^, Gi = l-7-ln(l + e)-^ii^^l-7: 



1 + e 

7 is the Euler's constant, and S{Q) is given by 

^zoo-f-A 



s{^) = r ^u^- = r ^v{e)e-''^(^^ , v{e) = J- 

J-ioo+A 27ri 27r smO 



(30) 
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The latter expression for S{^) in (|3l|) was obtained choosing the integration contour along the line 3u = in the u 
complex plane. Asymptotical expressions for S{^) may be easily obtained by the saddle-point method: 



s{n) 
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, as < 1. 
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Figure d shows nS {VL) as a function of In fJ. 
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The asymptotics of n(y) at small y is defined by the expressions After some simple calculations we have: 



vHv) 
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The two expressions, (|29|) and 
function in the integral in Eq. 



]), can be sewed together using the expression in the intermediate region, where the 
3) can be expanded up to the first order in a, ct- This yields in the region < y < 1: 
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(34) 



To establish regions of validity for three expression of the DF, let us set = 1 — A, A <C 1. Comparing Eq.(|34|) with 
Eqs.(|^ and (|3|), one can find, that Eq.(|^) is valid if A < Ai - a/eln (1/e), and Eq.(|^ — if A > A2 - eln (e/e). 
Let us note, that in spite of the low-conductivity cut-off for the universal distribution ( p9| ) is very close to the mean 
value in terms of the conductivity itself, Ai <C 1, which ensures central momenta of the conductivity to be of the 
order of a, this cut-off is small in terms of the universally fluctuating variable fJ, fJi ^ aP , p « 1 + 1/ In (1/e). 

The distribution function S (51) arise naturally in f d chain of random resistors, if to require scaling form of the 
distribution function of A-length chain specific resistivities: T (i?, A) = T{R\~°-), or Q (s, \) = Q (sA") in the Laplace 
representation. Then from Q{s,n\) — {s/n, X) one immediately has: Q (x) = exp [— Cx^/*^^+"'] . Evaluating 
its inverse Laplace T(r), and assuming a ^ 1, which is true in Id case, we have after the proper rescaling of the 
integration variable: 



a 



-l/a 



(35) 



which is essentially the same formula as Eqs. ( |29| , |30| ). 

It should be noted that, due to the nature of approximations used in the derivation of the RGMK equations, it 
is the distribution function at large resistivities (or low conductivities), which is most suspected to be inadeqately 
reproduced. Indeed, in reality some of the percolative paths at the threshold infinite cluster inevitably have return 
parts, where the diffusive particle moves in opposite direction. Such a paths are not taken into account in the RGMK 
scheme, which can be especially clearly demonstrated using hierarchical structures approach in the derivation 
This leads us e.g. to the overestimation of the conductivity value, resulting in lower value of the exponent a = t/v, 
which we found to be 1.890 vs. 2.25 ± 0.04 obtained through numerical simulation jl^ in 3d, and 0.818 vs. 0.98 in 
2d. Still the main part of conductivities distribution function, described by Eq. (p9|), seems to be valid in nearly- Id 
case. 

In conclusion let us note that the method suggested enables one to deal not only with distribution functions of 
conductivities in the percolative systems — it may also be applied to treat fluctuations of random variables in other 
disordered systems, e.g. ones described by random coupling Ising and Potts models. 
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FIG. 1. Limiting form of the distribution function for the rescaled logariphm of the conductivity. 
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